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L  PLASMA  THEORY  AND  SIMULATION 

A.  Lower  Hybrid  Drift  Instability;  2d  Simulations  and  Nonlocal  Theory. 

Yu-Jiuan  Chen  (Prof.  C.K.  Birdsall,  W.M.  Nevins,  LLNL) 

This  project  is  now  completed.  Two  ERL  reports  are  being  distributed  with  this  QPR. 
The  abstracts  are  given  in  Section  III-C. 

B.  Hybrid  Simulations  of  Theta  Pinch  Rotational  Instabilities. 

Douglas  Hamed  (Prof.  C.  K.  Birdsall) 

This  project  is  now  completed.  ERL  reports  are  being  distributed  with  this  QPR.  The 
abstracts  are  given  in  Section  III-C. 

C.  Stability  of  the  Thermal-Barrier  Cell  of  a  Tandem  Mirror:  Electrostatic  Simula¬ 
tions 


V.  Thomas  (Prof.  C.  K.  Birdsall) 

We  are  starting  to  investigate  the  stability  of  the  thermal-barrier  cell  of  a  tandem  mirror 
machine  to  electrostatic  ion-ion  two-stream  type  instabilities  using  a  1-D  electrostatic  quasineu¬ 
tral  particle  simulation  code.  The  code  model  uses  particle  ions  with  two  electron  components, 
a  "hot"  electron  component  and  a  "warm”  Boltzmann  component.  The  hot  electron  component 
is  modeled  as  a  fixed  negative  charge  density  confined  to  near  the  magnetic  field  minimum. 
Physically,  the  hot  electrons  represent  ECRH  heated  electrons  trapped  near  the  minimum  in 
the  magnetic  field  due  to  their  very  large  magnetic  moments.  The  warm  Boltzmann  electrons 
provide  charge  neutrality.  A  force  is  applied  to  the  ions  to  model  the  effects  of  the 

magnetic  well.  The  field  is  solved  for  in  this  code  using 

ne  T'  •  nt  —  nfh~j(B)  (1) 

Here  ne  T'  represents  the  warm  Boltzmann  electrons  and  n/ixtd(B)  represents  the  hot  elec¬ 
tron  component.  The  quantity  n,  is  accumulated  from  the  particle  ions. 

In  order  to  achieve  a  self-consistent  initialization  for  the  code  we  are  required  to  solve 
-».♦  -h* 

ne  r*  —  e  T‘  J  nl(vi,(i,^,B)dfidv:  -  n/^iB)  (2) 

for  4  as  a  function  of  r,  the  coordinate  along  the  field  line.  Then,  using  this  self-consistent 
potential,  it  is  possible  to  assign  appropriate  m  and  vi  values  to  the  ions. 

The  choice  assumed  for  the  ion  distribution  function  is  very  important  for  obtaining  well 
behaved  self-consistent  potentials*.  Some  distribution  functions  give  multi-valued  potentials 
near  the  mirror  throat  and  other  distribution  functions  generate  potentials  having  sheaths,  that 
is  locations  where  the  potential  has  a  discontinuity. 

We  are  using  a  distribution  function  given  in  the  aforementioned  reference  which  is 
based  upon  Fokker-Planck  studies  of  the  thermal-barrier  cells.  The  distribution  function  for 
ions  trapped  in  the  magnetic  and  electrostatic  well  of  the  thermal-barrier  cell  is  taken  to  be 


'Cohen,  R.  H.,  "Axial  Potential  Profiles  in  Thermal-Barrier  Cells,"  .Vuciaer  Ftaton  21.  219.  February  INI 
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f, 


Cexp 


E-anB^x 
(a-  or, 


and  the  distribution  function  for  the  passing  ions  is  given  in 


fp 


(3) 

(4) 


Here  C  is  a  normalization  constant,  E  is  the  energy  of  the  ion,  n  is  the  magnetic  moment  of 
the  ion  (which  is  treated  as  an  invariant),  and  is  the  magnetic  held  magnitude  at  the  mir¬ 
ror  throat.  The  quantity  a  is  a  parameter  which  determines  the  relative  densities  of  the  trapped 
and  passing  particles  in  the  Thermal- Barrier  cell.  An  important  feature  of  this  model  distribu¬ 
tion  is  that  the  distribution  function  is  continuous  across  the  separatrix.  This  leads  to  self- 
consistent  potentials  which  are,  in  general,  single  valued  and  which  vary  continuously  from  the 
mirror  throat  to  the  interior  of  the  thermal-barrier  ceil.  Distribution  functions  which  are  not 
continuous  across  the  separatrix  do  not  have  these  characteristics  in  geueral.  The  reader  is 
referred  to  the  reference  mentioned  previously  for  more  details. 


D.  Alfren  Ion-Cyclotron  Instability:  Physical  Description 


Niels  F.  Otani 


The  Alfven  ion-cyclotron  instability  (AIC)  is  a  velocity-space  instability  believed  to  occur 
in  the  ring-current  region  of  the  Earth's  magnetosphere  [Thorne,  1975]  and  in  certain  theta- 
pinch  [Davidson.  Ogden,  1975]  and  tandem  mirror  [Smith,  Nevins,  1981]  fusion  plasma  experi¬ 
ments.  As  is  generally  true  for  velocity-space  instabilities,  the  source  of  free  energy  responsible 
for  driving  the  instability  is  quite  evident;  in  this  case  it  is  Ta  >  T»  Somewhat  uncharacteris¬ 
tic  of  velocity-space  instabilities  however  is  the  fact  that  AIC  possesses  an  instability  mechan¬ 
ism  describabie  in  simple  physical  terms.  This  mechanism  was  described  briefly  in  Tqjima  and 
Dawson,  (1980]  and  will  be  amplified  here. 

We  note  that  the  assumptions  made  in  this  description  are  justified  by  the  consistency  of 
the  model.  In  particular  we  will  be  assuming  that  Tm  -  0  and  that  a*  =  tu„,  in  order  to  simplify 
the  description.  We  are  not  claiming  by  these  assumptions  that  AIC  does  not  occur  if  these 
assumptions  are  violated;  in  fact,  for  Tm  >  0,  AIC  is  strongest  when  <■>  is  appreciably  less  than 
[Davidson  and  Ogden,  1976].  It  will  be  clear  however  that  the  essence  of  the  mechanism  to 
be  described  holds  for  AIC  in  general. 

AIC  is  a  left-polarized,  electromagnetic  mode  which  typically  propagates  parallel  to  an 
unperturbed  magnetic  field  B0  which  we  assume  to  be  static  and  uniform.  We  also  assume 
<u  =  <uc/,  so  the  fields  and  electron  E|  x  Bo  drifts  (designated  (v^),)  in  Figure  1  are  rotating 
counterclockwise  around  the  axis  Bo  with  frequency  o>rt.  Since  the  local  unperturbed  ion  per¬ 
pendicular  velocity  distributions  (shown  as  disks)  are  also  rotating  with  their  natural  gyrofre- 
quency  <ocl,  we  may  think  of  Figure  1  as  rigidly  routing  with  frequency  tua,  apart  from  growth 
effects  due  to  the  insubility.  Alternatively,  we  may  regard  Figure  1  as  a  view  of  the  physics 
from  the  corouting  frame;  in  this  frame  everything  is  sutionary  except  that  the  perturbed 
quantities  exhibit  pure  growth  without  relative  change  in  orienution.  Note  however  that,  in 
this  view,  all  quantities  in  Figure  1  continue  to  be  lab  frame  quantities.  We  further  observe 
that  at  a  fixed  time  t,  a  translation  of  A/4  in  the  direction  of  propagation  produces  a  simple 
clockwise  roution  by  w/2  of  the  physical  picture. 

With  this  in  mind,  start  with  the  perturbed  magnetic  field  Bi  which,  as  shown  at  r  -  0  in 
Figure  1,  is  finite  since  this  is  an  electromagnetic  mode  and  perpendicular  to  Bo  since 
k  •  Bj  -  0.  The  quantity  dBi/dr  has  components  both  along  Bt  and  perpendicular  to  Bi  point¬ 
ing  in  the  direction  of  roution,  due  to  the  growth  rate  y  and  the  rigid  roution  w  =  oia,  respec¬ 
tively.  Faraday's  law  dicutes  that  the  wave  electric  field  E|  must  encircle  dBi/dr,  this  explains 
the  orienution  of  Ei  at  z  -  ±A/4  in  Figure  1. 
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This  electric  field  at  z  —  ±X/4  produces  two  contributions  to  the  perturbed  current.  First, 
the  electrons  drift  with  velocity  (yj)t  -  cE)XlV5o  in  the  direction  perpendicular  to  both  Et 
and  Bo,  assuming  negligible  electron  mass.  Second,  in  the  corotating  frame  of  a  given  ion,  the 
centrifugal  force  cancels  the  Lorentz  force  and  the  ion  is  simply  accelerated  in  the.  direction  of 
the  statically  oriented  electric  field.  This  is  pictorially  represented  in  Figure  1  at  r  -  ±A/4  by 
the  uncentered  positions  of  the  velocity  disk  relative  to  the  axis  of  rotation  Bo. 

A  third  contribution  to  the  current  at  z  -  ±  A/4  arises  from  the  interaction  of  Bt  with  the 
velocity  distribution  represented  in  Figure  1  at  z  -  0.  Ions  on  the  left  side  of  the  disk  are  mov¬ 
ing  to  the  left  and  consequently  feel  a  perturbing  Lorentz  force  ev0ix  Bt/  c  directed  antiparallel 
to  Bo-  Ions  on  the  right  side  of  the  disk  similarly  feel  a  force  parallel  to  Bo.  Since  the  orienta¬ 
tion  of  B]  is  fixed  relative  to  the  ions  instantaneously  on  the  left  (due  to  rigid  rotation)  they  are 
displaced  in  the  —Bo  direction.  Similarly,  rightward  moving  ions  are  displaced  in  the  +Bo  direc¬ 
tion.  These  displacements  are  represented  in  Figure  1  by  a  tilting  of  the  disk.  Thus  the  left¬ 
ward  moving  ions  are  preferentially  moved  toward  z  -  -A/4,  contributing  to  leftward  moving 
current.  We  see  this  contribution  to  the  current  tends  to  reinforce  Bi  at  z  -  0.  In  fact  this 
contribution  is  the  destabilizing  factor,  since  clearly  the  other  contributions  represent  rightward 
moving  current. 

In  summary,  three  contributions  to  the  perturbed  current  have  been  described.  The  effect 
of  two  of  these  contributions,  the  electron  and  ion  responses  to  the  perturbed  electric  field,  is 
stabilizing.  The  third  contribution,  resulting  from  the  selective  transport  of  ions  with  different 
perpendicular  velocities  up  and  down  the  zero-order  magnetic  held,  is  destabilizing  in  that  it 
tends  to  augment  the  perturbing  magnetic  field  which  caused  it. 
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E.  Relativistic  Quantum  Formulation  of  the  Boltzmann  Collision  Integral 

Kwang-Youl  Kim  (Professor  C.  K.  Birdsail) 

The  relativistic  quantum  conductivity  tensor  for  the  relativistic  plasma  is  derived  by  using 
the  relativistic  Liouville  equation.  The  quantum  relativistic  collision  integral  due  to  scattering 
via  covariant  photons  is  derived  by  using  relativistic  dielectric  response  function.  Our 
Boltzmann  collision  integral  includes  both  the  collective  fluctuation  effect  and  the  discreteness 
of  particles  on  the  same  footing.  The  corresponding  classical  results  are  obtained  by  taking  the 
appropriate  limits. 
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n.  CODE  DEVELOPMENT  AND  MAINTENANCE 

A.  Implicit  Particle  Simulation  Using  Moments  with  Orbit  Averaging 

.  V.  Thomas  (B.  I.  Cohen,  Prof.  C.  K.  Birdsail) 

The  project  is  now  completed.  A  report  has  been  written  (at  LLNL).  See  Section  III-D. 

B.  Plasma  Sheath  Simulation  in  Id 

S.  Rousset  (Prof.  C.  K.  Birdsail) 

The  numerical  model  described  in  the  previous  Quarterly  Progress  Report  has  been  imple¬ 
mented.  We  find  that  the  simulated  behaviour  of  the  sheath  agrees  in  some  respects  with 
well-known  phenomena  as  predicted  by  theory. 

A  new  subroutine  (called  ENTER)  has  been  added,  which  provides  for  the  simulation  of 
an  incoming  emitter  current.  In  its  present  version,  this  routine  adds  constant  numbers  of  par¬ 
ticles  at  each  time  step  (there  can  be  one  such  constant  number  per  species).  The  panicles 
enter  the  system  at  x— 0;  note  that  we  also  tried  to  let  new  panicles  enter  the  sheath  at  several 
different  positions,  but  this  did  not  seem  to  yield  any  measurable  improvement.  The  initial 
velocities  of  the  entering  particles  were  chosen  (almost)  randomly,  so  that  our  results  had  very 
high  noise  levels;  this  last  point  is  still  being  worked  upon. 

.  Since  we  have  new  particles  entering  the  system  at  each  time  step,  and  since  some  anays 
in  the  code  have  sizes  proportional  to  the  total  number  of  simulated  panicles,  we  had  to  deal 
with  the  problem  of  defining  the  sizes  of  the  arrays  in  a  way  such  that  we  never  attempt  to  load 
panicles  beyond  the  limits  of  these  arrays.  On  the  one  hand,  a  straightforward  solution  is  to 
repack  the  arrays  from  time  to  time,  and  to  expand  them  whenever  we  try  to  go  above  their 
upper  bounds.  On  the  other  hand,  we  could  declare  the  sizes  of  the  arrays  to  be  much  larger 
than  the  initial  number  of  particles,  adding  new  panicles  at  the  end  of  the  area  (of  the  anay) 
already  occupied  by  other  panicles.  The  first  solution  is  very  costly,  and  the  second  one  sets  a 
very  stringent  limit  on  the  number  of  time  steps  in  one  run.  We  chose  the  following  comprom¬ 
ise:  the  size  of  the  arrays  are  chosen  equal  to  the  maximum  number  of  panicles  that  will  ever 
be  simultaneously  in  the  system;  the  arrays  are  repacked  at  each  time  step.  This  very  last  point 
means  that  we  could  not  keep  the  vectorized  option  in  subroutine  MOVE,  but  we  believe  that 
our  solution  is  most  efficient  in  terms  of  total  bank  time. 

A  few  other  technical  features  have  been  implemented:  the  plots  do  not  include  the 
inverted  part  of  the  sheath  any  more;  a  call  to  the  PRIOR  subroutine  has  been  added,  thus 
allowing  the  priority  to  be  set  automatically  at  run  time. 

As  a  consequence  of  the  modifications  described  above  and  in  the  previous  QPR,  several 
new  parameters  can  now  be  given  in  the  input  file: 


Name  Function  Default 

For  the  whole  system 

hiprio  Maximum  bid  priority  allowed  for  the  job  1.0 

I  Physical  length  of  the  system  (excluding  inverted  sheath)  tr 

ng  Number  of  spatial  cells  in  the  true  sheath  only  16 

nprie  If  nprto— 0.  then  disable  the  PRIOR  routine  0 

nrank  Position  of  the  job  in  the  priority  queue  20 

phiO  Amplitude  of  the  externally  applied  potential  0.0 

For  each  species 

a  Number  of  particles  loaded  at  time  zero  64 

nenter  Number  of  new  panicles  loaded  at  each  time  step  1 

amax  Maximum  number  of  particles  simultaneously  in  the  system  64 


We  found  the  new  version  of  the  code  to  reproduce  some  of  the  physical  behaviours 
predicted  for  a  sheath.  In  particular,  some  features  of  Child's  law  have  been  successfully 
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checked  for  three  regimes:  without  space  charge;  with  space  charge  and  without  initial  velocity 
for  particles  at  the  emitter;  and  with  space  charge  and  with  initial  velocity. 

References 
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J.  K.  Boyd,  "PRIOR,  a  Subroutine  to  Efficiently  Set  the  Running  Priority  on  the  Cray  Computer,"  Memorandum 
number  MFE/CP/80-2432t,  1980. 

C.  Quiet-Start  Method  Comparisons: 

Random,  Bit-Reversed,  and  Fibonacci  Numbers 

Yu-Jiuan  Chen,  and  H.  Stephen  Au-Yeung 
(Prof.  C.  K.  Birdsall;  Dr.  A.  B.  Langdon,  LLNL) 

In  QPR  IV,  1980,  we  showed  that  the  Fibonacci  number  and  the  bit-reversed  methods 
produce  very  good  uniformity  and  that  the  random  number  generator  produces  considerable 
local  bunching.  As  a  further  test,  we  have  also  used  the  cumulative  velocity  distribution  func¬ 
tion  to  assign  velocities  vx  to  a  set  of  beams  whose  envelope  is  /(vx)  (as  in  ESI)  and  then  used 
the  three  different  number  sets  to  scramble  the  particle  positions,  x,  to  reduce  correlations. 

For  a  given  non-drifting  Maxwellian  velocity  distribution  function,  /( vx)  —  e  1  ' ,  the 

’goodness”  of  the  local  thermal  equilibria  was  examined,  by  calculating  the  average  of  the  first 
three  velocity  moments  over  various  widths  8x  <  L  (.L  is  the  system  length  and  v,  is  the  ther¬ 
mal  speed  of  this  velocity  distribution  function).  By  varying  the  particle  number  V,  we  then 
compared  the  width  5x  at  which  the  calculated  average  local  velocity  moments  departed  by  10% 
from  the  expected  values,  for  each  of  the  three  different  number  sets,  as  shown  in  Figure  1. 
Note  that  the  expected  odd  velocity  moments  vanish  for  a  non-drifting  Maxwellian  velocity  dis¬ 
tribution.  Therefore,  fix’s  were  measured  at  which  the  average  local  odd  velocity  moments 
were  10%  of  v,"  where  m  —  1  in  Figure  1(a)  and  3  in  Figure  1(c),  respectively.  We  found  that 
for  a  large  particle  number  N  (^  103)  both  the  Fibonacci  number  method  and  bit-reversed 
method  produce  very  good  local  thermal  equilibria. 


Figure  1.  Goodness  of  the  local  thermal  equilibrium  velocity  moments  as  a  function  of  the  number  of  simulation  parti¬ 
cles  used  .V  For  random,  bit-reversed  and  Fibonacci  methods,  ix  were  measured  at  which  (a)  <  v>»,/  v,  -  10% 
(b)  K  vJ>,I  -  v*>/  .,J  -  10%  and  (c)  <»*>,,/  v,J  -  10%.  Here  <  »*>„  is  the  calculated  average  m*  velocity 
moment  over  width  ix. 
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Journal  publication:  (Completed  while  Friedman  was  here;  distributed  with  last  QPR) 

A.  Friedman,  R.N.  Sudan,  and  J.  Denavit,  "A  Linearized  3D  Hybrid  Code  for  Stability  Studies 
of  Field-Reversed  Ion  Ring,"  J.  Comput.  Phys.  40,  1-35,  March  1981. 
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A  Linearized  3D  Hybrid  Code  for  Stability  Studies 
of  Field-Reversed  Ion  Rings 

A.  Friedman*  and  R.  N.  Sudan 

Laboratory  of  Plasma  Studies,  Cornell  University,  Ithaca,  New  York  14853 

AND 

J.  Denavit 

Department  of  Mechanical  Engineering  and  Astronomical  Sciences. 
Northwestern  University,  Evanston,  Illinois  60201 

Received  Miy  22,  1979;  reviled  February  7,  1980 


A  linearized  3D  hybrid  code,  R1NGHYBRID,  for  the  study  of  the  low-frequency  stability 
of  high-beta  equilibria  it  described.  In  particular,  this  code  is  suitable  for  application  to  field- 
reversed  ion  rings  and  mirror  plasmas.  The  code  models  the  development  in  time  of  pertur¬ 
bations  having  a  specified  azimuthal  mode  number,  I,  about  an  axisymmctric  equilibrium 
state.  The  equilibrium  current  is  purely  azimuthal  and  is  carried  entirely  by  the  energetic  ion 
component,  which  is  modeled  by  discrete  simulation  particles.  This  current  is  consistent  with 
the  equilibrium  magnetic  field  and  hot-component  charge  density.  In  addition  to  the  hot-ion 
component,  a  cold  uniform  density  background  of  ions  and  a  complement  of  massless 
electrons  are  described  by  fluid  equations  which  incorporate  a  scalar  resistivity.  Also 
presented  are  verifications  of  code  performance,  including  tests  of  cold-plasma  normal  modes 
and  of  plasma  return  currents  arising  in  response  to  an  azimuthal  “drive'’  current  rising 
linearly  with  time. 


I.  Introduction 

In  recent  years  there  has  been  a  rebirth  of  interest  in  field-reversed  configurations, 
which  have  the  advantage  of  possessing  a  region  of  closed  lines  for  high  0  plasma 
confinement  surrounded  by  a  region  of  open  lines.  This  topology  requires  a  relatively 
simple  design  for  the  external  magnetic  field  coils  that  leads  to  a  major  simplification 
in  fusion)  reactor  design.  The  region  of  closed  lines  is  produced  by  internal  currents  in 
the  plasma.  When  such  currents  are  force-free  the  device  is  known  as  a  spheromak 
( 1 1.  If  the  plasma  currents  flow  in  the  azimuthal  direction,  only  poloidal  fields  are 
present  and  we  have  a  reversed  field  9-pinch  [2J.  A  reversed-field  mirror  (3)  is  a 
configuration  similar  to  the  9-pinch  except  for  a  lack  of  axisymmetry  introduced  by 
having  minimum  B  quadrupole  type  coils  to  provide  the  external  field.  The  build-up 
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It  is  found  that  a  finite  perturbation  of  the  ion  orbits  lends  to  a  nonlinear  frequency  shift  that  reduces  the 
mode  frequency  and  has  a  weak  stabilizing  effect  on  the  lower-hybrid  drift  instability.  This  result  is  obtained 
from  a  self-consistent  solution  of  the  Vlaaov-Poisson  equations  using  perturbation  theory  in  which  the 
nonlinear  dielectric  function  and  the  nonlinear  temporal  evolution  of  a  single  unstable  mode  in  the  low  drift 
velocity  regime  are  calculated  analytically.  The  nonlinear  frequency  shift  does  not  seem  to  be  a  potent 
saturation  mechanism  in  a  collisionless  plasma,  but  may  be  more  relevant  when  there  are  ion-ion  collisions. 


I.  INTRODUCTION 

The  linear  theory  of  the  lower-hybrid  drift  instability 
is  well  understood  and  has  been  discussed  in  detail  by 
Davidson  el  al.1  When  the  amplitude  of  the  wave  is 
small  but  finite  after  a  time  equal  to  many  multiples  of  the 
growth  time,  further  evolution  will  be  different  from 
exponential  growth  at  very  small  amplitudes.  In  order 
to  analyze  this,  the  nonlinear  dielectric  response  func¬ 
tion  and  the  nonlinear  tempora1  evolution  of  a  single 
unstable  mode  are  derived  self-consistently  by  using 
perturbation  theory  to  solve  the  Vlasov  and  Poisson 
equations.  The  single-mode  approximation  is  valid  for 
the  instability  when  the  plasma  parameters  are  close 
to  those  for  linear  marginal  stability.  This  requires 
vE  «  vti  (the  low  drift  velocity  regime)  for  the  lower 
hybrid  drift  instability,  where  ve*cE/B  is  the  equilib¬ 
rium  E  x  B  drift  velocity  and  uf,=  T/M  is  the  ion  ther¬ 
mal  speed.  A/  and  Tt  are  the  Ion  mass  and  tempera¬ 
ture,  respectively.  Similar  single  mode  studies  in¬ 
clude;  modulation  of  the  Langmuir  wave  due  to  weak 
nonlinearity2  nonlinear  evolution  of  drift- cyclotron 
and  drift-cone  instabilities  both  In  theory4'1  and  simu¬ 
lation.' 


In  this  paper  we  demonstrate  that  a  finite  amplitude 
perturbation  of  the  ion  orbits  develops  during  growth 
which  leads  to  a  weakly  stabilizing  nonlinear  shift  of 
the  lower-hybrid  drift  mode  frequency.  The  largest 
shifts  in  the  mode  frequency  occur  for  modes  with 
wavelengths  much  longer  than  that  of  the  most  unstable 
mode  (h  at  which  wavenumbers  the  lower-hybrid 

drift  instability  converts  into  the  drift  cyclotron  insta¬ 
bility.  For  simplicity,  we  use  a  one -dimensional  slab 
configuration  shown  in  Fig.  1;  wave  propagation  is  in 
the  x  direction;  the  magnetic  field  is  uniform  and  in  the 
z  direction;  the  density  gradient  is  in  the  y  direction. 
The  ions  are  treated  as  unmagnetized  because  the  wave 
frequency  and  growth  rate  are  much  greater  than  the 
ion  cyclotron  frequency.  The  ions  are  in  force  balance; 
that  is ,  the  force  due  to  the  ion  pressure  gradient  in  they 
di  rection  cancels  that  of  the  equilibrium  electric  field. 

The  lower-hybrid  drift  instability  is  analyzed  in  the  elec¬ 
trostatic  limit;  electromagnetic  effects  are  assumed 
to  be  small. 

In  Sec.  n,  the  nonlinear  dielectric  response  is  calcu¬ 
lated  by  solving  the  coupled  Vlasov-Poisson  equations, 
and  a  nonlinear  dispersion  relation  is  obtained.  Sec¬ 
tion  III  is  devoted  to  a  derivation  of  the  time  evolution 
of  the  lower -hybrid  drift  instability.  The  field  energy 
level  at  which  saturation  might  occur  and  the  frequency 
shift  due  to  the  finite  amplitude  of  the  wave  are  also 
determined.  Ftnaliy,  conclusions  and  a  comparative 
discussion  of  several  saturation  mechanisms  are  given 
in  Sec.  IV. 
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B. 

Abstracts  follow  of  papers  submitted  to  APS  Division  of  Plasma  Physics  meeting,  October 
12-16,  1981,  New  York 


Nonlocal  Properties  of  the  Lower  Hybrid  Drift 
Instabilities*  YU- J 113 AN  CHEN,  UC  Berkeley,  W.M.  NEVINS 
LLNL  Using  Kinetic  theory  in  the  low  drift  velocity 
regime  and  considering  only  electrostatic  perturbations 
in  slab  geometry,  we  found  that  the  lower  hybrid  drift 
instability  is  a  negative  energy  wave  driven  by  reso¬ 
nant  ions  wherever  sgn(ui/ky)=sgn  v*e(ky)  as  predicted 
by  the  local  theory,  which  also  says  that  the  local 
growth  rate  is  largest  in  the  region  where  the  rel¬ 
ative  drift  velocity  of  electrons  and  ions  is  largest. 
In  our  non-local  analysis  this  region  is  an  energy 
source  for  normal  modes.  The  wave  energy  propagates 
from  this  region  across  the  zero  order  density  gradient 
to  regions  where  the  electron  drift  velocity  vE  equals 
the  wave  phase  velocity  (u/ky)  and  is  then  absorbed 
by  the  resonant  electrons.  As  a  result,  the  unstable 
normal  modes  are  localized  between  the  regions  of 
high  relative  drift  and  electron  resonance.  The  exis¬ 
tence  of  propagation  along  Vn  was  observed  first  in 
2-D  particle  simulations  motivating  this  theory;  the 
two  now  check  very  well.  These  simulations  show  that 
the  modes  saturate  due  to  both  EbcB  trapping  and  cur¬ 
rent  relaxation.  *Research  partially  supported  by  ONR 
Contract  //N00014-77-C-05 78  (Berkeley) ,  and  in  part  by 
DOE  Contract  #  W-7405-ENG-48  (Livermore). 

Rotational  Instabilities  in  the  Field-Reversed  Theta  Pinch  — 
Results  of  Hybrid  Simulations  ,  D.S.  HARNED,  U.C.  Berke¬ 
ley.*  A  2d  quasineutra!  hybrid  simulation  code  has  been  used 
to  study  rotational  instabilities  in  the  field-reversed  theta  pinch. 
Ions  are  particles  so  that  both  the  stabilizing  effects  of  finite 
Larmor  radius  and  the  destabilizing  effects  of  resonant  ions  are 
included.  For  the  non-reversed  case,  our  results  for  the  m— 2 
instability,  where  m  is  the  azimuthal  mode  number,  are  com¬ 
parable  to  previous  theory1,2.  In  contrast,  for  the  reversed 
case,  using  parameters  like  those  of  FRX-B,  we  have  observed 
instability  for  a<l  where  Sle  and  fl,  are 

the  electron  and  ion  rotational  frequencies).  This  new  result  is 
consistent  with  observations  in  FRX-B.  The  threshold  value 
has  been  found  to  depend  critically  on  the  shape  of  the  pres¬ 
sure  profile.  As  the  profile  becomes  more  peaked,  the  thres¬ 
hold  decreases.  For  sufficiently  peaked  profiles,  unstable 
modes  with  m>  2  have  been  observed.  Nonlinearity  has  been 
found  to  reduce  the  growth  rates  for  the  m— 2  instability  at 
large  amplitude. 

1.  J.P.  Freidberg  and  L.D.  Pearlstein,  Phys.  Fluids  21,  1207 
(1978). 

2.  C.E.  Seyler,  Phys.  Fluids  22,  2324  (1979). 

*  Work  supported  by  ONR  Contract  No.  N00014-77-C0578. 
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C. 

Abstracts  follow  of  recent  ERL  reports,  most  intended  for  submission  to  journals,  as 
noted;  distributed  with  this  QPR.  (Parts  of  Chen  and  Harned  Ph.D.  dissertations) 


REVIEW  OF  THE  LOWER  HYBRID  DRIFT  INSTABILITY 
AND  ITS  SATURATION  MECHANISMS 

Yu-Jiuan  Chen 

Electronics  Research  Laboratory 
University  of  California,  Berkeley,  CA.  94720 

ABSTRACT 

A  survey  is  given  of  previous  research  on  the  linear  theory  of  the  lower 
hybrid  drift  instability  (Sec.  I).  A  summary  is  made  of  research  recently  com¬ 
pleted  as  well  as  research  currently  in  progress  on  the  nonlinear  theory  and 
computer  simulations  of  the  lower  hybrid  drift  instability  (Sec.  11). 


Memorandum  No.  UCB/ERL  M81/60  20  August  1981 
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STABILIZATION  OF  THE  LOWER  HYBRID  DRIFT 
INSTABILITY  BY  RESONANT  ELECTRONS 

Yu-Jiuan  Chen 

Electronics  Research  Laboratory 
University  of  California,  Berkeley,  California  94720 

William  M.  Nevins 

Lawrence  Livermore  National  Laboratory 
University  of  California,  Livermore,  California  94550 

Charles  K.  Birdsall 

Electronics  Research  Laboratory 
University  of  California,  Berkeley,  California  94  720 

ABSTRACT 

The  lower  hybrid  drift  instability  was  studied  with  a  two  dimensional  elec¬ 
trostatic  simulation  code.  Simulations  showed  good  agreement  of  the  measured 
local  growth  rates  and  frequencies  with  the  results  of  local  theory  during  the 
early  stage  of  wave  growth.  At  later  times  nonlocal  effects  become  important, 
and  a  coherent  mode  structure  develops.  This  normal  mode  was  observed  to 
propagate  up  the  density  gradient. 

At  zero  plasma  beta  and  zero  electron  temperature,  we  found  that  the 
lower  hybrid  drift  instability  is  stabilized  by  the  local  current  relaxation  due  to 
both  ion  quasilinear  diffusion  and  electron  £x  B  trapping  which  causes  electron 
heating  to  occur. 

ERL  Memo.  No.  M81/72  16  September  1981. 

Prepared  for  submittal  to  Physics  of  Fluids. 
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Quasineutral  Hybrid  Simulation  of  Macroscopic  Plasma 

Phenomena 

Douglas  S.  Homed 

Electronics  Research  Laboratory 
University  of  California 
Berkeley,  CA  94720 

ABSTRACT 

A  method  for  solving  the  quasineutral  hybrid  plasma  equations  in 
two  dimensions  is  presented,  using  full  ion  dynamics  and  inertialess 
electrons.  The  method  is  extended  to  allow  plasma-vacuum  interfaces 
of  arbitrary  shape.  The  algorithm  is  applied  to  the  study  of  rotational 
instabilities  in  theta  pinch  Vlasov  equilibria. 

ERL  Memo.  No.  M81/71  16  September  1981. 

Prepared  for  submittal  to  Journal  of  Computational  Physics. 


Kink  Instabilities  in  Long  Ion  Layers 

Douglas  S.  Homed 

Electronics  Research  Laboratory 
University  of  California 
Berkeley,  CA  94720 

ABSTRACT 

Kink  instabilities  in  long  ion  layers  immersed  in  a  dense  back* 
ground  plasma  are  studied.  A  numerical  extension  of  the  analytic 
model  of  Lovelace  indicates  that  these  instabilities  will  occur  for  values 
of  the  self-magnetic  field  index  below  those  predicted  previously.  A 
quasineutral  hybrid  simulation  code  has  been  used  to  verify  these  lower 
thresholds.  The  simulations  also  show  that  the  end  of  exponential 
growth  occurs  due  to  a  nonlinear  shift  in  the  betatron  frequency  at  large 
amplitude,  producing  an  increase  in  layer  thickness  and  a  layer  which 
has  many  non-axis-encircling  ions. 

ERL  Memo.  No.  M81/73  16  September  1981. 

Prepared  for  submittal  to  Physics  of  Fluids. 


Rotational  Instabilities  in  the  Field-Reversed  Theta  Pinch: 


Results  of  Hybrid  Simulations 


Douglas  S.  Harried 

Electronics  Research  Laboratory 
University  of  California 
Berkeley,  C A  94720 

ABSTRACT 

Rotational  instabilities  in  rigidly  rotating  field -reversed  theta  pinch 
equilibria  are  studied  using  a  quasineutral  hybrid  simulation  code.  We 
observe  unstable  m-2  modes  at  levels  of  ion  rotation  below  instability 
thresholds  predicted  by  finite  Larmor  radius  fluid  theory.  Nonlinear 
effects  are  found  to  reduce  the  growth  rate  and  lower  the  real  frequency 
of  the  m— 2  mode  at  large  amplitude.  Instabilities  with  m>  2  have  been 
observed  for  some  strongly  reversed  cases.  It  is  also  found  that  growth 
rates  for  these  instabilities  can  be  greatly  reduced  by  increasing  the  ratio 
of  the  plasma  radius  to  the  ion  Larmor  radius. 

ERL  Memo.  No.  M81/70.  16  September  1981. 

Prepared  for  submittal  to  Nuclear  Fusion. 
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D. 

Abstract  of  Report  UCRL-86429  at  LLNL  follows.  This  was  prepared  for  submission  to 
Journal  of  Computational  Physics  (work  done  at  LLNL  and  here). 

B.I.  Cohen,  R.P.  Freis,  V.  Thomas,  "Orbit-Averaged  Implicit  Particle  Codes" 

ORBIT-AVERAGED  IMPLICIT  PARTICLE  CODES 

Bruce  I.  Cohen  and  Robert  P.  Freis 
Lawrence  Livermore  National  Laboratory,  University  of  California 
Livermore,  California  94550 

and 

Vincent  Thomas 

Electronics  Research  Laboratory,  University  of  California 
Berkeley,  California  94720 

ABSTRACT 

This  paper  reports  the  merging  of  orbit-averaged  particle  code 
techniques  with  recently  developed  implicit  methods  to  perform  numerically 
stable  and  accurate  particle  simulations.  Implicitness  relaxes  time-step 
constraints.  Orbit  averaging  allows  one  ensemble  of  particles  to  represent 
many  ensembles  of  particles  distributed  with  respect  to  a  phase  variable  over 
the  phase-space  trajectories.  In  simulation  models  where  the  electromagnetic 
fields  are  calculated  with  the  time  step  much  longer  than  the  step  used  for 
particle  pushing,  a  large  reduction  in  the  statistical  requirements  on 
particles  can  result'  from  orbit  averaging.  Implicitness  and  orbit  averaging 
can  extend  the  applicability  of  particle  codes  to  the  simulation  of  long 
time-scale  plasma  physics  phenomena  such  as  transport  and  low-frequency 
instabilities  in  controlled  fusion  devices.  Difference  equations  for 
electrostatic  and  magneto-inductive  physics  models  are  presented,  and  analyses 
of  the  numerical  stability  of  each  scheme  are  given.  Simulation  examples  are 
presented  for  a  one-dimensional  electrostatic  model  and  for  a  two-dimensional 
magneto-inductive  model. 

UCRL  -  86429  PREPRINT  (LLNL)  29  July  1981 
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